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Abstract. We show that the inner regions of a cool accretion disk in an X-ray binary can transform into an 
advective, ion supported accretion flow (an optically thin ADAF, here called ISAF), through events involving 
only the known properties of the Coulomb interaction in a two-temperature plasma, standard radiation processes, 
and viscous heating. The optically thin inner edge of the disk is heated to a few 100 keV by the strong flux of hot 
ions from the surrounding hot ISAF. We show that he resident ions in this 'warm' disk are thermally unstable due 
to internal viscous heating, and heat up to their virial temperature. The innermost disk regions thus evaporate 
and feed the ISAF. These processes are demonstrated with time dependent calculations of a two-temperature 
plasma in vertical hydrostatic equilibrium, including heating by external ions, internal proton-electron energy 
exchange, and viscous heating. The process complements the 'coronal' evaporation mechanism which operates at 
larger distances from the central object. 
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1. Introduction 

Accreting galactic black holes (BHC) and active galactic 
nuclei (AGN) are often observed with two different spec- 
tral components: a soft component which is probably due 
to a multi-color blackbody from an optically thick, geo- 



metrically thin standard disk (Shakura & Sunyaev 1973, 
SSD), and a hard component which is linked to an op- 
tically thin and geometrically thick flow. The hard com- 
ponent (an approximate power law with high energy cut 
off at E c w 100 keV) is most likely produced by in- 
verse Co mpton scattering of soft ph ot ons on a hot therma l 
plasma ( Sunyaev fc Titarchuk 198C ). Shapiro et al. (1976 ) 
showed that accretion can take place in the form of a two- 
temperature plasma, with the properties needed to pro- 
duce such Comptonized radiation. Stabl e, optically thin , 
two-tcmpcraturc flows were studied by Ichimaru (1977) 
and Rees et al. (1982). Extensive theoretical work on these 
accretion flows, stressing the role of the advection of in- 
ternal energy was done by|Narayan fc Yi (1994 |l995a| H), 
for a review see Narayan et al. (1998), see also Esin et al. 
(1997) and references therein. 

The effects of advection are the same for the optically 
thin two-temperature flows and optically thick, radiation 
supported flows, and both types of flow are now customar- 
ily called 'ADAF's. The distinction between these types, 
which was explicit in the older labels 'ion supported' and 
'radiation supported', is still needed in many applications, 
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however. For this reason, we propose here the name 'ion 
supported accretion flow' or ISAF as a means of distin- 
guishing the optically thin type of ADAF from the radia- 
tion supported ('RSAF') type. 

ISAFs have attracted attention because of their poten- 
tial to explain the spectra of X-ray transients (Esin et al. 
1997 and references therein). The observations are con- 
sistent with an accretion flow that consists of two zones: 
an interior ISAF that extends from the black hole hori- 
zon to a transition radius i?tr, followed by an optically 
thick, geometrically thin and cool standard disk outside 
i?tr- A partial overlap between the two regions is probable 
since observations show evidence for the close vicinity of 
hot and cold matter in the central regions of BHCs and 
AGNs. This is indicated by a K Q iron fluorescence line 
at 6.4 keV and a Compton reflection component between 
« 10-30 keV. 

A critical element of such an accretion geometry is 
the change from the geometrically thin SSD to the hot 
ISAF flow at the transition radius. An alternative inter- 
pretation of the spectra is given by the 'magnetic flare' or 
disk corona model (Haardt & Maraschi 1991 , Maraschi 
& Haardt 1997 and references therein, di Matteo et al. 
1999, Merloni et al. 2000 and references therein), in which 
this transition is absent and the X-ray emitting plasma 
is a layer on top of a disk, heated by mechanical energy 
transfer from the disk. 

Signatures of a disk transition radius have recently 
been found in the power density spectra of Cyg X- 



2 



H.C. Spruit and B. Deufel: The transition from a cool disk to an ion supported flow 




Fig. 1. Energy channels in a cool accretion disk embedded in a hot region. The energy dissipation in the cool disk is 
assumed small compared the atmosphere (a corona or an ion supported ADAF, here called ISAF). Squared boxes show 
physical processes, round boxes the particles involved. Heavy arrows and boxes show the main energy channel: viscous 
dissipation in the ISAF heats the protons there, which illuminate the cool disk below. By Coulomb interactions the 
ISAF protons loose their energy mainly to electrons, producing a warm (~ 80 keV) layer which radiates this energy 
by Compton-upscattering of soft photons from the cool disk below. The protons in the warm surface layer are largely 
outside this main energy channel. Near the transition radius R m mtr the cool part disappears and the warm part heats 
up to several hundred keV. Then the energy exchange of disk protons with the electrons is slow (shown by the saw 
tooth line) and viscous heating of protons becomes important. 



l QChurazov ct al. 200l| ), ( |Gilfanov fc Revnivtsev 2000]) 
From an analysis of the reflection component Di Salvo 



et al. (2001) showed that in Cyg X-l the transition ra- 
dius R tI is located between 10i?s and 70i?s, if the ob- 
served reflected spectrum is due to a smeared component, 
or 6i?s < -Rtr < 20i?s if the reflection is unsmeared, e.g. 
from the companion star or the outer disk. 

We take these observations as reasonable indications 
that a transition from a cool optically thick to a hot op- 
tically thin accretion flow does in fact take place. But 
how the SSD-ISAF transition works is still under debate 



(cf M anmoto et al. 2000 and references therein). [Meyer 



et al. [2000J) propose that the transition from the cold 



disk to the optically thin flow is due to a heat flow by 
electron conduction from a hot, friction-heated corona to 



the cold disk below (see also Meyer fe Meyer-Hofmcister 
1994[)7]rhis model has its maximum evaporation efficiency 
at a large distance from the hole (a few 100 Schwarzschild 
radii). If, due to a high accretion rate in the cool disk, 
not all material is evaporated until that distance, the cool 
disk will survive until the last stable orbit. A transition 
radius further in than 100 Schwarzschild radii is inconsis- 
tent with this picture. Rozanska fe Czcrny (2000| ) inves- 
tigate conductive and radiative coupling of an accretion 
powered corona with an underlying cool disk. For low ac- 
cretion rates they find that the disk completely evapo- 
rates whereas high accretion rates prevent the SSD-ISAF 
transition as in Meyer et al. (2000|) . From a mathematical 



point of view Abramowicz et al. (1998| ) show that, if the 



transition region is not too wide, the region must rotate 
with super Kcplerian orbital speed. Based on this prop- 
erty Kato fe Manmoto (2000| ) demonstrate that trapped 
low-frequency oscillations are possible in the transition 
region. 



In this paper we show how the inner disk regions, where 
the coronal evaporation process does not work, can evap- 
orate into an ISAF. Only a few, well-known ingredients 
need to be invoked: the coupling of protons and electrons 
by the Coulomb interaction in a fully ionized plasma, stan- 
dard radiation processes, and viscous heating. 

The starting point is the view, supported by the ob- 
servations mentioned, that an ISAF and a cool disk can 
coexist. That is, there is a partial overlap between the cold 
disk and the ISAF. In the overlap region, there is a very 
strong interaction between the two, since the ISAF con- 
sists of energetic ions (10-100 MeV) that penetrate the 
cool disk to a significant depth. The goal in the next sec- 
tions is to determine the nature of the energy and mass 
exchange in this interaction region, and to show that it 
will lead to evaporation at the inner edge of the cool disk. 
The argument is then closed by determining the conditions 
under which these processes can consistently lead to the 
coexistence of the disk and the ISAF that was assumed at 
the outset. 



1.1. Interaction between ISAF and cool disk 

The essence of the processes described below is given by 
the energetic interaction between ISAF and disk. In Fig. 
|l] we have sketched the main energy channels involved in 
this interaction. The ultimate source of energy is the re- 
lease of gravitational binding energy. We assume here that 
a significant fraction (at least a few tens of per cent) of 
this energy goes into the protons (on account of their much 
higher mass compared to the electrons). This viscous en- 
ergy release predominantly takes place in the hot region 
(the ISAF). The protons and electrons in the ISAF form 
a two-temperature plasma, where the protons have tern- 
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peratures near their virial temperature. The electrons are 
much cooler due to their strong interaction with the radia- 
tion field and the slow rate at which they exchange energy 
with the ISAF protons. 

Due to their low temperature, the conductive energy 
flux carried from the ISAF to the disk by the electrons is 
negligible. This is one of the reasons why the coronal evap- 
oration mechanism that functions well at larger distances 
from the central mass fails for a two-temperature plasma. 
The energy flux to the cool disk is carried almost entirely 
by the ions. This ion energy flux can not take over the role 
played by electron conduction in the coronal evaporation 
process, for two reasons. One is that the mean free path 
of the ions is not negligible as it is in the case of electrons. 
The stopping length of 10-100 MeV ions penetrating into 
a cool disk corresponds to a Thomson optical depth of 
order unity. This prevents the development of the thin en- 
ergy deposition layer that is needed to heat the plasma 
to high temperatures. More important is the fact that the 
energy carried by these ions primarily heats the electrons 
of the disk. These in turn radiate it very effectively by 
inverse Compton scattering and limit the plasma temper- 
ature in the interaction region to values of ~ 100 keV, well 
below the virial temperature. Instead of evaporation, the 
loss of the ISAF ions to the cool disk is a very effective 
condensation process. This condensation is an important 
sink to an ISAF flow generated by evaporation at larger 
distances, and has to be overcome by a sufficiently power- 
ful evaporation process at some location in the disk. This 
location will turn out to be a region near the inner edge 
of the disk. 



1.2. A 'warm' surface layer on the cool disk 

The ISAF and the cool disk are separated by a 'warm' sur- 
face skin of temperature ~ 100 keV, intermediate between 
the ion temperature of the hot ISAF and that of the cool 
disk. This heated surface layer is produced by the energy 
flux from the ISAF to the cool disk. This flux can be in the 
form of radiation (photons) o r particles (ions ). Radiative 



coup ling was investigated by Haardt fc Maraschi (1991 
1993| ) in the context of their 'two-phase model' for the 



hard X-ray spectra of accreting black holes. Explanation 
of these spectra as resulting from the coupling between an 
ISAF and a cool disk via hot the protons was proposed by 
Spruit (1997| ) and jBpruit fc Haardt (2000| ), and studied in 



greater detail by Deufel fc Spruit (200C| ) [henceforth pa- 
per I ] and Deufel ct al. (2001a| ) [henceforth paper II ] . The 
physics of this process has been studied before in the con- 



text of accretion onto a neutron star surface by Zel'dovich 
fc Shakura (1969| ); |Alme fc Wilson (1973| ), pcufcl ct al 
(2001b|) . 



The surface layer radiates its energy by Compton- 
upscattering of soft photons from the cool disk below. The 
energy of the ISAF is thus carried to the disk by the hot 
ISAF protons, and radiated away by the disk electrons. 



The protons in the warm surface layer are largely outside 
this main energy channel. 

In the present analysis we concentrate on ion heating, 
which is necessarily very strong in the overlap region be- 
tween an ion supported flow and a cool disk. A radiative 
flux from the ISAF can exist additionally, but in order to 
separate these contributions we assume here that the ISAF 
itself is radiativcly inefficient, so that its direct contribu- 
tion to the radiation from the accretion flow negligible. 
Finally, the viscous energy release in the warm layer turns 
out to be small compared with incident energy flux from 
the ISAF. Viscous heating is thus energetically unimpor- 
tant in the surface layer, but it will turn out to be crucial 
to the evaporation process. 

Finally there is the cool disk (~ 1 keV) below the warm 
layer. This region is outside of the reach of the impinging 
hot protons. The cool disk serves as a thermalizer for the 
downward directed radiation of the warm layer. It is the 
source of the soft photons which are Comptonized in the 
warm region and keep it at moderate temperatures. The 
viscous energy release in the cool disk is assumed to be 
negligible as a source of soft photons. 

In paper II we have shown that close to the inner edge 
of an accretion disk, where the surface density (and the 
optical depth) of the disk gets small, the penetration of 
virialized protons heats the entire vertical disk structure 
to temperatures of several 100 keV (which is equivalent 
to a disappearance of the cold part in Fig. |l|). This corre- 
sponds to the first step in our model and is at the same 
time the starting point of the present investigation. At 
high temperatures the time scale for establishing thermal 
equilibrium between the disk protons and electrons is not 
short compared to the dynamical time scale any more. The 
viscous energy channel (due to internal heating of protons 
by friction) gets important now because coupling to the 
electrons is weak. 

It will turn out that the physics described depends 
almost only on the dimensionless radius from the hole 
r = B,/Rs, the dimensionless accretion rate M/M^dd, and 
the dimensionless viscosity parameter a. The results are 
thus scalcable between AGN and BHC cases. Where ex- 
plicit values of the physical parameters are needed, we 
take those of a typical BHC case. 

1.3. Coulomb interactions in an ionized plasma 

In the next three subsections, we briefly review the physics 
associated with the penetration of the protons into the 
disk. This has been discussed before in detail in the refer- 
ences given in the previous subsection. Here, we address a 
few conceptual issues, such as the validity of the approx- 
imation that the energy of the incident protons is trans- 
ferred mainly to the electrons in the disk, the accuracy of a 
nonrelativistic treatment, and the charge balance between 
ISAF and disk. 

At the typical energies of the protons incident on the 
cool disk, the energy loss is mostly by long-range Coulomb 
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interactions with the electrons in the disk (small-angle 
scattering on the large number of electrons in a Debye 
sphere). This is opposite to the case of protons with a 
temperature near that of the plasma in which they move. 
In the latter case, the equilibration among the protons 
is faster than between electrons and protons, by a factor 
of order (m-p/nic) 1 / 2 . To see how this apparent contra- 
diction is resolved, consider the basic result for the en- 
ergy loss of a charged particle moving in a fully ionized, 
charge-neutral plasma. This was derived by Spitzer (1962) 
(making use of Chandrasekhar's (1942) earlier results on 
dynamical friction) . Introduce as a measure of distance in 
the plasma the Thomson optical depth r, i.e. dr = aTndl, 
where or = 8ire / (3m 2 c ) is the Thomson cross section, n 
the electron density and / the distance. The rate of change 
of energy E = m p v 2 /2 of a proton moving with velocity 
v in a field of particles with charge e and mass rrif (the 
'field particles' in Spitzer's nomenclature) is then given by 
Spitzer's Eq. 5-15. In our notation, this can be written in 
terms of the energy loss length Tf for interaction with the 
field particles f , 

1 ,dE s 



In this limit, the energy loss is thus predominantly to the 
electrons. A related case is that of the ionization losses of 
fast particles in neutral matter (for references see Ryter et 
al. 1970). The case of an ionized plasma is simpler, since 
the electrons are not bound in atoms. The change from 
proton-dominated loss to electron-dominated loss takes 
place at an intermediate velocity v c , at which F(x p ) ~ 1 
but x e still small, so that F(x c ) « Equating r e and 

r p then yields 



(6) 



(7) 



For the electron temperatures we encounter in our models, 
T ~ 100 keV, energy loss to the field protons can thus be 
neglected for incoming protons with energy E ^ 1 MeV. 
This is the case in all calculations presented here. 



(— )f = 31nA(— ) 2 (-) 4 (l H ^-)F(xf), (1) 1.4. Corrections at high and low energies 



E'dr 



rrit 



where In A is the Coulomb logarithm (which is determined 
by the size of the Debye sphere) . Here 



F(x) = i/j(x) — xip'(x) 7 



(2) 



where tj) is the error function and ip' its derivative, and 
X{ = v[mi/ {2kT)] 1 / 2 is (up to a numerical factor) the ratio 
of the incident proton's velocity to the thermal velocity of 
the field particles. The limiting forms of F are 



F(x) -> x 3 /3 (x 0), 



F 



1 (a; — * oo). 



(3) 



We can evaluate ([!]) under the assumption that the field 
particles that are most relevant for the energy loss are the 
protons or the electrons, respectively, and compare the loss 
lengths. If the incident proton has velocity comparable 
with the thermal velocity of the field protons, we have 
x p ~ 1, and x c — (r7i c /m p ) 1 / 2 a:p <C 1. F(x p ) is then of 
order unity and F(x c ) « x^/S. Setting f — e respectively 



f = p in (|lj) and taking the ratio, we have 



(4) 



The loss length for interaction with the electrons is thus 
much longer than for interaction with the protons, and 
the interaction with electrons can be neglected. This is 
the well known result for the relaxation of a proton dis- 
tribution in a plasma that is not too far from its thermal 
equilibrium. 

For incoming protons of high energy, however, the re- 
sult is different because x c is not sufficiently small any 
more. In the high-z; limit, F(x e ) — F(x p ) 
has 



1, and one 



— » 2— = 10 



-3 



(5) 



Spitzer's treatment is non-relativistic, while virialized 
ISAF protons near the hole can reach sub-relativistic tem- 
peratures. A fully relativistic treatment of the Coulomb 
interactions in a plasma has been given by Stepney & 
Guilbert (1983). We have compared the classical treat- 
ment according to Spitzer's theory with this relativistic 
result in Deufel et al. (2001), and found it to be accurate 
to better than 5% for proton temperatures < 100 MeV. 
The classical approximation in Spitzer's analysis therefore 
does not introduce a significant error for the problem con- 
sidered here. 

For high energies, the Coulomb energy loss becomes so 
small that loss by direct nuclear collisions becomes com- 
petitive. This happens (cf. Stepney & Guilbert 1983) at 
E <; 300 MeV, an energy that is not reached by virialized 
protons except in the tail of their distribution. We ignore 
these direct nuclear collisions. Note, however, that a grad- 
ual nuclear processing by such collisions can be important 
(Aharonian & Sunyaev, 1984), in particular for the pro- 
duction of the Lithium. The Lithium overabundances seen 
in the companions of LMXB (Martin at al., 1994a), may 
in fact be a characteristic signature of the interaction of 
an ADAF and a disk described here (Martin et al. 1994b, 
Spruit 1997). 

As the protons slow down, they eventually equilibrate 
with the field protons. This last part of the process is not 
accurately described by the energy loss formula (|TJ). In 
addition to the simple energy loss of a particle moving on 
a straight path through the plasma, one has to take into 
account the random drift in direction and energy resulting 
from the interaction with the fluctuating electric field in 
the plasma. This drift can be ignored to first order (end 
of Sect. 5.2 in Spitzer 1962), but takes over in the final 
process of equilibration with the plasma. This last phase 
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involves negligible energy transfer compared with the ini- 
tial energy of the protons in the present calculations, and 
can be ignored here. 

1.5. Charge balance 

The protons penetrating into the disk imply a current that 
has to be balanced by a 'return current'. As in all such 
situations, this return current results from the electric field 
that builds up due to the proton current. As this electric 
field develops, it drives a flow of electrons from the ADAF 
to the disk which maintains the charge balance. Since the 
electron density in the disk is high, the return current does 
not involve a high field strength. 



2. The evaporation process 

Any mechanism heating the protons in the disk, even if 



The energy exchange timescale by Coulomb interac- 
tions between nonrelativistic protons and electrons with 
Maxwellian distributions at temperatures T p and T c , re- 
spectively, is given by Spitzer's (1962) classical result: 



t cp = — — ; , js " P . . ^ 253 ' 



8(2ttto c ) 1 /2 e 4 n e In A 



In A 



s, 



(8) 



where m p , m are the mass of the proton and electron, 
respectively, k is the Boltzmann constant, e the charge 
of the electron, n c the electron density and In A « 20 
the Coulomb logarithm. The approximate equality is used 
since we consider cases where initially the electron tem- 
perature is of the order of the proton temperature and 
therefore the contribution of can safely be neglected. 

The relativistic correction to this result is small for the 



only small, can potentially give rise to a thermal mstabil- 
ity since the disk protons do not loose energy efficiently. 
Their energy loss by radiation is negligible on account of 
their high mass, and the transfer of energy to the disk elec- 
trons by Coulomb interaction is inefficient in a sufficiently 
hot plasma. In addition, by the nature of the Coulomb 
interaction, the time scale for the disk protons to equili- 
brate wit h th e disk electrons increases with temperature 
(cf Sect. F3). Depending on the strength of the mecha- 
nism heating the protons and its temperature dependence, 
a runaway may occur in which the disk protons continue 
to heat up as their cooling rate continues to decrease. 

One heating process to be taken into account is inter- 
nal viscous energy release (due to friction). At the low sur- 
face densities in the region of interest, the amount of en- 
ergy released by friction is small compared with the heat- 
ing flux by the external ISAF protons. If we may assume 
that a reasonable fraction of the viscous energy release 
goes into the protons, this small amount can still be quite 
important for the energy balance of the disk protons be- 
cause of their low energy loss rate, and must therefore 
be taken into account. A second heating process is en- 
ergy transfer from the incident ISAF protons to protons 
in the disk by Coulomb interaction. Though we have seen 



conditions of interest here. In Deufel et al. (2001a ) we have 
shown that Spitzer's formalism is accurate to better than 
5% compared to the relativistic treatment by Stepney & 



Guilbert (1983), for proton energies below 100 MeV and 



electron temperatures kT c < 50 keV. 

Next, the viscous heating timescale in an accretion disk 
is given by 



tth 



1 



(9) 



where fl = [GMj i? 3 ) 1 / 2 is the Kepler angular velocity and 
a is the viscosity parameter (shakura & Sunyaev 1973). 
Comparing Eqs. flpl) yields a critical electron density 



3m p fc 3/2 



8\/27r m c e 4 In A 



■T? /2 afl . 



(10) 



( Sect. 1.3) that the incident protons loose their energy 
mainly to the electrons, the small amount transferred to 
the disk protons may conceivably be relevant for their en- 
ergy balance. In the following we investigate these heating 
mechanisms in turn, an conclude that viscous heating can 
give rise to a runaway, while heating of the disk protons by 
Coulomb interaction with the incident protons does not. 

2.1. Viscous heating of the disk protons, instability 

The energy balance of the protons in the warm disk is con- 
veniently described in terms of heating and cooling time 
scales due to the processes involved. First we consider the 
cooling time scale by transfer of energy from protons to 
electrons in the disk. 



For electron densities n c < n* viscous heating in an ac- 
cretion disk works faster than proton-electron coupling. 
In those regions the protons can not loose their energy 
fast enough, and heat up. This increases n*. At the same 
time, the higher proton temperature causes the layer to 
expand, by hydrostatic equilibrium. On both accounts, 
the Coulomb coupling between protons and electrons de- 
creases, and the heating of the protons accelerates. Thus 
we expect viscous heating of the protons to lead to ther- 
mal instability wherever the density is less than given by 
Eq. (|l(]). Since the density decreases steeply with height 
in the atmosphere, there is always a level above which it 
is unstable. A new equilibrium is reached only when the 
protons reach the virial temperature. The unstable part 
of the atmosphere has then expanded to H/R ~ 1, and it 
is then, effectively, part of the ISAF in which the disk is 
embedded. The unstable part of the atmosphere has evap- 
orated, feeding the ISAF. Whether such an instability has 
an effect on the global accretion properties depends on the 
mass in the unstable region. 

The amount of mass involved in the instability is the 
mass in the atmosphere above the level where the density 
has dropped to n%. Let this mass be M* (per unit surface 
area of the disk) . The time scale on which the protons heat 
up in the unstable region is just the viscous heating time 
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scale t t h- The rate at which the atmosphere evaporates is 
thus approximately 



M « atlM* 



(11) 



The unstable mass M* is found by integrating the den- 
sity in the atmosphere upward from the level z* where 



/>oo 

M* = / m p n c (z) dz — m p n* c H\/2f(i 



(12) 



where H is the scale height of the disk atmosphere 
(assumed isothermal) with density stratification, n c = 
n exp[— u 2 ] with u 2 = z 2 /2H 2 . The function f(u*) is 



/(«*) 



exp [u 



*2 



I du = e % 
1 2 



(l-erf[w*]) (13) 



To evaluate Eq. (|T|) we need to know the dimension- 
less critical height above the midplane, u* . The level at 
which n c ~ n* depends on the surface density of the disk, 
being higher at large surface density. We measure the sur- 
face density by the Thomson scattering optical depth T X j 2l 
measured from the midplane to the surface of the disk. In 
terms of the density at the midplane n , this is: 



n/2 = n (JtH \J ir/2 . 

With Eq. © we can then calculate u*: 



1/2 ,2 



-In 



;t m p ' k 



aT 2 



V2m 5 c /2 c 4 In A T i/2 



1/2 



(14) 



(15) 



Above a certain temperature (for a given a) u* is not 
defined, as the value within the square brackets in Eq. ( |l5|) 
drops below zero. This means that the whole atmosphere 
is subject to the instability, and we can set u* = 0. In this 
cases f{u*) = \/tt/2 (the maximum value for /). 

Using Eq. (|l2|) to replace M* in Eq. ([n|) this yields, 
with Eq. @: 



3/2 



M 



27T raj 2 e 4 In A 



a 2 riT 2 f(u*) 



(16) 



sj2 k T/nip is the isothermal 



Here H = c s /f2 and c 
sound speed. 

An evaporating part exists in every disk atmosphere, 
but at low temperatures only the highest layers evaporate. 
The mass loss is unimportant in such a case and can not 
change the properties of the disk. Eq. (|l^) shows that the 
evaporation rate is a relatively strong function of R, a and 
T c . Whether situations exist where this evaporation rate 
becomes important is investigated in Sect. 



2.2. Heating by proton-proton interactions 

In this subsection we show that Coulomb interaction of 
the incident protons with the resident disk protons is rel- 
atively unimportant for the conditions encountered. It is 
included in the numerical results reported in the next sec- 
tion, however. 

The main energy channel is from the hot protons to the 



electrons in the warm surface layer (see Sect. 1.3 and pa- 
per II) . The properties of the Coulomb interactions at the 
temperature of the incoming protons are such that only a 
small fraction of their energy is transferred to the protons 
in the layer (which we will call 'field protons' here follow- 
ing the terminology in Spitzcr 1962). The energy budget 
of the field protons is also small, however (cf. Fig. |l|), so 
the effect of this channel on the proton temperature in the 
warm layer has to be considered. 

The interaction of a fast proton with the protons in a 
much cooler plasma is given by Spitzer's result discussed 
in Sect. |l|. In terms of the rate of change of the energy of 
the incoming proton, it can be expressed as 



d£y 
dt 



87re 4 



n p In A [ip{x) — xip'(x)] , 



(17) 



where n p is the density of the (cool) field protons. The 
subscript p' indicates the hot penetrating protons and p 
the field protons. Here ip( x ) an( i V''( a; ) are the error func- 
tion and its derivative and x 2 = m p v 2 /2kTis the ratio of 
the velocity of the incoming protons to the thermal veloc- 
ity of the field protons. In the following we set ip — xip' = 1. 
This approximation is valid if the velocity of the incom- 
ing protons much exceeds the velocity of the field protons. 
In our case we consider virialized protons penetrating into 
the considerably cooler disk plasma, so the approximation 
is valid. 

The heating rate W p ' p per unit volume due to the in- 
teraction of hot protons with the field protons is then 



^p'p 



dE p , 
dt 



8^e 4 

7T in A n„ n„> v p > 

m p vp 



(18) 



This can be compared with the rate of energy transfer 
from the field protons to the electrons in the warm layer 



as given by Eq. (5-30) from Spitzer (1962) 



Wp C = 



'-Cp 



8(2tt to c ) 1/2 e 4 n 2 In A kT c - kT p 
3^ (fcT c )3/ 2 



cp- 



(19) 



The 



where we have used Eq. (|8|) to replace t, 
approximate equality indicates that the contribution 
from m c /m p T p can again be neglected. Using v p > = 
yj2kT p * /m p , we compare Eqs. ( |l8| , ^9|) and obtain an esti- 
mate for the equilibrium temperature of the disk protons 
exposed to penetrating hot protons: 
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T n « T P + 114 



n r . 



(20) 



If the mass flux in the ISAF and its viscosity parameter 
a are of the same order as in the cool disk, a thin disk 
approximation for both flows yields n p </n c ~ (T e /T p /) 3 / 2 . 
In the inner region of the accretion flow, where T c ~ 100 
keV and T p > > lOMeV, the second term in @) is then 
negligible. 

We conclude that heating of the ions in the warm disk 
by the incident ions alone is not sufficient to drive the ion 
temperature away from equilibrium with the electrons. For 
this to happen, viscous heating has to be included. 

2.3. Numerical simulation of the evaporation 

We test the inferred heating instability numerically, start- 
ing with a warm disk of moderate optical depth and in- 
vestigating its temporal evolution due to internal viscous 
heating, and p'-p and p-e energy exchange as described 
in the previous subsection. 

We use a plane-parallel, one-dimensional model geom- 
etry. The vertical density distribution through the atmo- 
sphere is found from the equation of hydrostatic equilib- 
rium as in paper II, except that here we do not account 
for the pressure exerted by the decelerating protons. 



Th 



temporal evolution oi tte lemperalure profile is 
computed with a fourth order Runge-KulU niolliod. The 

timestep At of each integration is set to the shortest 
timescale of the heating and cooling processes involved. 
The change of temperature per timestep within a certain 
volume is due to the rate of change of the enthalpy there, 
A t w = pc p A t T. 

The temperature increase per timestep as a function 
of optical depth due to viscous heating of the protons is 
then given by 



Ar vi sc(r) 9 1 
-^— = -an e (r)kT p (r)n K .-^- 



(21) 



where c p is the specific heat at constant pressure. 

The temperature change of the field protons due to the 
impinging hot protons is computed similar to the proton 
illumination method described in detail in paper I and II. 
Thus we follow protons from a Maxwellian distribution 
with virial temperature through the warm atmosphere of 
the disk and record their energy losses. The energy flux in 
the illuminating protons is related to the mass condensa- 
tion rate [see Eq. ( p6[ ) and Sect. |]ior details]. The kinetic 
energy decrease of the fast protons (p') due to their inter- 
actions with the field protons (p) as a function of optical 
depth can be derived from Eq. (|l^). This gives the heat- 
ing rate A p / p and the contribution of p'p exchange to the 
temperature change per time step is 



AIV p (t) _ A 



p'p 



At 



(22) 



Fig. 2. Evolution of the temperature profile due to the vis- 
cous instability of the protons. The model layer has optical 
depth T1/2 = 1.5 and initial temperature To iP = Xb.c = 500 
keV (dashed line). Dashed-dotted lines show the profile 
after 10 thermal time scales, solid lines after t=20. Values 
of the viscosity parameter a are 0.1, 0.3 and 0.5 (numbers 
at the lines). For comparison the temperature profile af- 
ter t = 20 £th due to p'p interactions alone is also shown 
(dotted line). 

Finally the rate of change of the temperature due to 
the electron-proton [ep] coupling is [according to Bpitzer 



(1962|) 1 



AT ep (r) _ T„(t) - T p (r) 



At 



t e p(r) 



(23) 



To complete the model the electron temperature in the 
disk needs to be specified, which is regulated by the radia- 
tion processes. We have considered these in some detail in 
paper II, where we found that for the present conditions 
electron temperatures of 100-500 keV result. The radia- 
tion losses by inverse Compton scattering increase quite 
rapidly with T c and with the increasing optical depth when 
pair production sets in. This makes the electron temper- 
ature relatively insensitive. For the present purpose, it is 
sufficient to keep the electrons at a constant temperature, 
which we treat as a parameter of the model. 

When the disk protons in our simulation have reached 
their local virial temperature, we do not allow a further 
temperature increase. In a more realistic calculation, this 
limit on the temperature would come about through the 
advection of internal energy. Such a calculation requires 
a more detailed model of the ISAF and its sources and 
sinks of mass and energy, which is beyond the scope of 
the present calculation. 

As an example we consider a warm disk with an optical 
depth (from the surface to the midplane) T-yi^ = 1.5, an 
initially isothermal temperature profile T c = T p = 500 
keV, illuminated by an ISAF with accretion rate M\ = 
O.lMEdd, where -MEdd = 4wcR+m p /<JTh is the Eddington 
accretion rate. 
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Fig. 3. Local evaporation and conden- 
sation rates from a warm (100-500keV) 
layer of optical depth T1/2 = 1-5, 
formed by proton illumination of an ini- 
tially cool disk. Solid line shows the 
condensation rate onto the layer from 
an ISAF accreting at O.IMe, as a func- 
tion of distance from the hole. Dashed 
lines show the evaporation rates from 
the layer, for three values of its tem- 
perature. In the innermost region, the 
net effect is evaporation, feeding the 
ISAF, further out the ISAF condenses 
onto the layer. At lower ISAF accretion 
rates ( O.OIA/e, dotted) the evaporat- 
ing region becomes wider (see text for 
details) . 



Fig. ^| shows the temperature profiles for various val- 
ues of the viscosity parameter a after t = 10 t t h (dashed- 
dotted lines) and t = 20t t h (solid lines), where i t h = 
l/(af2) is the thermal time scale of the disk. For a = 0.5 
all disk protons heat up to the virial temperature, whereas 
for a = 0.3 and a = 0.1 only a part of the disk is sub- 
ject to the instability. The extent of the unstable par t is 
in excellent agreement with our estimate in Sect. 2.1. In 
test calculations in which only p'p heating and ep energy 
coupling arc included, the temperature increase of the disk 
protons is considerably smaller than with the viscous heat- 
ing included. This confirms our conclusion from Sect. [2.2| 
that p'p interactions alone can not heat the bulk of the 
internal protons considerably above the ambient electron 
temperature [cf. Eq. (pfj|)]. 

Summarizing our investigation of the energy channels 
in a warm disk we conclude that viscous energy release in 
a hot plasma causes a runaway temperature increase in 
the unstable upper part of the disk atmosphere, until the 
protons have reached their local virial temperature. The 
process affects a significant part of the stratification if the 
electron temperature is above ~ 100 keV. 

3. Evaporation vs. condensation rates 

In Sect. H we have shown that the internal viscous heat- 
ing of the disk protons leads to a mass evaporation rate 
according to Eq. (|l6|). 

At this point we do not yet know whether the mass loss 
Ai from a warm disk region is high enough to completely 
evaporate the disk. At the same time as the upper atmo- 
sphere of the warm disk evaporates, the hot protons from 
ISAF condense into it and increase the surface density. 
For an effective evaporation of the disk the mass loss rate 



must be higher than the condensation rate. To compute 
the condensation rate we need an estimate of the density 
in the ISAF. 

In our previous numerical simulations of warm disks 
(paper II) we parameterized proton mass flux from the 
ISAF by scaling the energy flux of the incident protons 
with the local energy dissipation rate in the ISAF. Here 
we adopt a slightly more realistic mass flux rate in an 
ISAF. 

In a thin disk approximation the surface density £1 of 
the ISAF, with accretion rate Mj, is 



3irac 2 



Mi[l-(3R S /R) 



1/21 



(24) 



We assume that the protons in the ISAF have a 
Maxwellian velocity distribution according to their local 
virial temperature. The mass flux rate </> [in gem -2 sec -1 ] 
from a Maxwellian through a surface is given by 



= Pp 



1 kT p 
2tt m r . 



(25) 



From Eq. (J2J) we can calculate the proton density p\ — 
Y>i/(2H\) in the ISAF. This yields the local mass flux rate 
onto the cool disk for protons at their virial temperature, 
i.e. the condensation rate: 



Mi 



91 = — 



8tT 3 / 2 R 2 



[1-{3R S /R) 



1/2] 



(26) 



We can now compare the mass condensation rate with 
the evaporation rate according to Eq. ([l6|). The evapora- 
tion rate is oc (aT) 2 . Therefore high disk temperatures 
and high a's favor mass evaporation. In paper II we have 
calculated these temperatures with a radiative transfer 
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Fig. 4. Decrease of the proton mass flux with depth in a 
thin, hydrostatic, isothermal disk with total optical depth 
Ttot = 2 ■ T1/2 = 3, at disk temperatures of 1, 50 and 200 
keV, for R = 8Rs and a = 0.1. At temperatures above 100 
keV, only a fraction of the incident protons is captured by 
the disk 



calculation that includes brems photon production, 77 
pair production and Comptonization. We found there that 
the temperature at the inner edge of a cool accretion disk 
reaches several 100 keV ('warm state'), while its optical 
depth T1/2 (at the border with the cool part of the disk) is 
of the order unity for accretion rates around O.lMEdd- We 
assume here that at a certain radius the warm disk state 
exists and, ignoring its radial structure, compare the local 
evaporation and condensation rates. 

Fig. |^ shows this comparison for different values for the 
viscosity parameter a, different temperatures and accre- 
tion rates. For values of the viscosity parameter a > 0.3 
and M = O.lMEdd evaporation dominates when the tem- 
perature of the warm state exceeds ~ 300 keV. As the 
ISAF density (i.e. the accretion rate) decreases, the con- 
densation rate of the protons into the the disk decreases 
and evaporation dominates over the condensation rates 
over a wider range of radii. 

The condensation rate given by Eq. ( p(f ) is actually an 
overestimate, since it assumes that all incident protons are 
stopped in the disk. While this is correct for cool disks, 
for a high temperature plasma the rate of the electron- 
proton energy exchange is small and a disk with low sur- 
face density gets optically thin for the penetrating hot 
protons. This is demonstrated by Fig. (5) where we show 
how the incident proton flux changes with depth into a 
warm disk. The Thomson optical depth in this example is 
(ftot = 2 • T1/2 — 3). At the temperatures of a cool stan- 
dard disk (< 1 keV) almost all protons are absorbed in this 
layer. But at the high temperatures of the warm state the 
disk is optically thin and practically all ISAF protons fly 
through the disk without being absorbed. The penetrat- 
ing protons do not add to the surface density in this case, 



and evaporation should therefore be possible even at lower 
values of a. 

At temperatures below 100 keV, the evaporation rate 
from the warm state into the ISAF is quite small and 
can not balance the loss from it by condensation. This 
is roughly the temperature of the warm surface layer on 
a cool optically thick disk heated by proton illumination 
(the situation sketched in Fig 1). Thus the relative impor- 
tance of evaporation and condensation reverses just at the 
point where the cool component disappears. As long as a 
cool disk is present, the thermal instability in its warm sur- 
face layer is relatively weak, while it effectively absorbs all 
incoming protons. Once the cool component is gone, the 
temperature and evaporation rate goes up, while at the 
same time mass condensation by stopping of protons in 
the disk becomes ineffective. 

3.1. Dependence on the accretion rate 

Figs. ||,|] show that net evaporation takes place close to 
the hole preferentially at low accretion rates in the ISAF, 
which may sound counterintuitive. It is a consequence of 
the fact that the evaporation rate does not depend on 
the flux of incident protons (cf. Eq. [If]). This again is a 
result of the fact that the field protons in the warm disk 
lie 'outside the main energy channel': the incident proton 
energy flux sets the electron temperature (through the 
Comptonization balance) but does not affect the proton 
temperature directly. The temperature of the field protons 
is determined by the secondary balance between viscous 
heating of the field protons and their energy loss to the 
electrons. 

At low accretion rates, however, the hot proton flux 
eventually becomes insufficient to keep the layer 'warm': 
it will cool down to low (~ 1 keV) temperatures by 
bremsstrahlung losses. We have studied this transition in 
paper II, where we found that a layer of optical depth 
7"i/2 = 1.5 can just be kept in the warm state for an ac- 
cretion rate in the ISAF of about 0.1 Me- 

At lower accretion rates, a warm disk is possible only 
when the bremsstrahlung losses are lower, at lower elec- 
tron densities. Since the electron density is proportional to 
the optical depth r of the disk, and the brems losses pro- 
portional to n 2 , the minimum ISAF accretion rate needed 
to maintain a warm disk state by proton illumination 
scales as r 2 . The optical depth of the disk vanishes to- 
wards its inner edge, so we expect that there is always a 
region close to the inner edge of the disk where evapora- 
tion takes place, even at very low accretion rates. 

Summarizing this argument, the optical depth r m of 
a warm disk at the point where it matches onto the cool 
disk, depends on the accretion rate in the ISAF as 

T -'-^" 3 {^) ■ (27) 
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4. An evaporating cool disk 

4.1. Thin disks with evaporation/condensation 

In the above we have started with the assumption that 
an ISAF and a cool disk coexist, and found the condi- 
tions under which the disk can feed mass into the ISAF. 
We have done this by considering the conditions at each 
distance from the hole separately. To turn the ingredients 
into a consistent picture, we have to consider the mass 
flux through the system, so that conditions as a func- 
tion of distances from the hole are connected, in the way 
sketched in Fig. 0. At large distance, we have a standard 
cool disk. Closer in, an ISAF surrounds it and condenses 
onto it, producing what we have called here the warm, 
proton-illuminated layer. Even further in, the vertical op- 
tical depth of the disk becomes too low to sustain cooling 
by brems losses, and the whole disk transforms to a 'warm 
disk' state (paper II) of nearly uniform temperature. The 
upper layers of this warm disk evaporate to feed the ISAF 
assumed at the outset, and at the inner edge of the disk, at 
some radius the entire mass flux through the disk has 
evaporated into the ISAF. We now investigate what the 
conditions are for such a radial structure to be possible in 
a steady state. 

The disk, whether in a cool (< 1 keV) or warm (~ 300 
keV) state, is still very cold compared to the local virial 
temperature, so that the standard thin disk approximation 
is valid. The difference with a standard steady disk is that 
the mass flux is now a function of distance, due to the 
condensation and evaporation from/to the ISAF. We first 
consider the modifications to the a-disk diffusion equation 
that this causes. 

The surface mass density of the disk material is 



E = 



pdz w 2po H , 



(28) 



where poi H denote the density at the midplane and the 
scale height of the cool disk, respectively. 

The change of the mass accretion rate with radius 
through the cool disk due to evaporation can be expressed 

by 



dM 

7m 



AttRM , 



(29) 



where we include an additional factor 2 to account for the 
mass loss on both sides of the disk. Integrating this equa- 
tion yields the mass loss (in g sec -1 ) due to evaporation 
in a ring between R\ and R, 



M fiv = 4tt / RMdR 



(30) 



The conservation of mass in a cool disk including the 
evaporation term [Eq. (|29|)] can be expressed by 



-(iffi) + QxiRXvR) + RM = , 



(31) 



where Ur is the radial drift velocity of the disk material. 

As in the standard derivation, the thin disk diffusion 
equation follows from the angular momentum equation, 
which now includes a term for the angular momentum 
carried with the evaporating/condensing material. For the 
present purpose, it is sufficient to assume that condensing 
and evaporating material just has the same specific angu- 
lar momentum, QkR 2 = as the disk. 

The equation for the angular momentum balance in an 
evaporating disk then is 



^-(REVlR 2 ) + -^-(KSvrQR 2 ) + RM SIR 2 
ot oR 

dR {SR dR ] ' 

where 

poo 

S = pvdz w YjV . 



(32) 



(33) 



We use the usual a prescription for the viscosity 
(ghakura fc Sunyaev 1975), 



v = a 



ft ■ 



(34) 



The approximate equality in Eq. fl33| ) holds if v can be 
considered independent of the geometrical height z. 

4.2. Steady state 

Multiplying the continuity equation [Eq. j3f|)1 with R 2 Q 
yields after subtraction from Eq. ( |32| ) an expression for 
the mass flux M in the cool disk in the stationary case 
(d/dt = 0): 



M = 6irR 1 / 2 -^(vXR 1 / 2 ) . 



(35) 



The thin disk equations, Eqs. ( |3l| , |32[) can be used in 
their time dependent form. For example, the evaporating 
inner regions could expand inward towards the last sta- 
ble orbit, or outwards. As a first application, however, we 
are interested in steady-state conditions. In the station- 
ary case we can use Eq. (|35| ) to replace M in Eq. (f29|). 
Integrating this expression yields 



R 1/2 J^(v^R 1/2 ) = i / RMdR + Ci 



dR 



(36) 



We integrate one more time and obtain the surface 
density distribution: 



2 f R dR 



Co 



RMdR + 2d + -^- , (37) 



3i?V2 J R; R i/2 R l/2 



where Ci, C2 are integration constants. Let the total mass 
loss due to evaporation from i?; to infinity be M eVtOC (cf. 
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Fig. 5. Evaporation from the warm disk: a cool disk (black, ~ 1 keV) partly extends into an ISAF and is exposed to 
the hot virialized protons from the hot torus (open arrows, light grey). Above and below the cool disk a heated surface 
layer is produced (~ 80 keV, dark grey) due to proton illumination. At the inner edge at Rq the surface density of 
the cool disk gets small and a warm disk develops (~ 300 — 500 keV, grey). Net mass loss into the ISAF evaporates 
the warm region between Ri < R < Rq (filled arrows). Due to the high temperatures in the warm disk this region is 
a source for hard X-rays. 



eq.|g). With Eq 

Ci, 



Ci = -L(M 

D7T 



Mo 







we find the integration constant 



(38) 



As in the standard treatment, the second integration 
constant is fixed by considering the conditions at the in- 
ner edge of the disk. This could be the last stable orbit if 
the disk extends all the way down to the hole, but more 
interesting is the case when the entire mass flux has evap- 
orated into an ISAF before reaching the hole. Thus we now 
assume that a steady state exists with the inner edge at 
some as yet unspecified distance The requirement of a 
steady state will then impose a condition on the parame- 
ters of the system that has to be satisfied. This will turn 
out to be a condition on the accretion rate. At the end 
of the calculation, we thus obtain a relation between the 
mass accretion rate and the position of the inner edge. 

Thus we set uHi = at the inner edge of the disk, 
as in standard accretion theory. With Eq. ( |30| , [37|j3S| ) and 
R = Ri we get an expression for the second integration 
constant C2, 



Co 



R 



1/2 



3vr 



■(M — M e v,oo) ■ 



(39) 



With the previous results the surface density distri- 
bution as a function of radius in a disk with evaporation 
losses is now 



if the position of the inner edge changes slowly. We are in- 
terested in a true stationary case, however. In this case the 
accretion rate in the disk will be exactly equal to the total 
mass loss due to evaporation in the disk, M — M 0ViOO , 
or in other words all matter drifting inward through the 
cool disk has evaporated when Ri is reached. Eq. (36) then 
simplifies to 



1 



dR 



torRWJn. R 1 / 



M ev (R). 



(41) 



4.3. Radial extent of the warm disk 



We can now estimate the distance over which the pro- 
cess of evaporation into the ISAF takes place. Let Ro be 
the innermost radius where the cool disk component ex- 
ists. Inside this, there is only a warm disk (see Fig. |J). 
Evaporation takes place both from the warm layer on top 
of the cool disk at R > Ro and from the warm disk re- 
gion Ri < R < Rq, but the warm disk region is expected 
to contribute most, since its temperature is significantly 
higher. Thus we equate, for the present approximate pur- 
pose, the evaporating region with the warm disk region. 
Assume that the relative extent 60 = Ro/Ri — 1 of the 
warm disk is small. The evaporation rate Eq. ([l^), with 
O Qi, depends only on temperature. The temperature 
of the warm disk is relatively fixed (paper II), so we can 
set M. constant as well. Eq. jj]] ) for the surface density 
as a function of distance from the inner edge is then 



1 



67^/2 



dR 



M ev (R) 



(40) 



3^ M 



[l-(Ri/R) 



1/2] 



From this expression the classical formula for the sur- 
face density in a cool disk can be recovered if one sets 
M cv = 0, i.e. when no evaporation takes place. 

This expression for the surface density, though strictly 
derived for steady conditions, is still approximately valid 



dR 



Ri 



£ = -A S 2 , 



where 



A 



■MR? 



MdRtz-M (R-Rif 



Ri 



(42) 



(43) 



(44) 
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In the warm state, little net condensation takes place 
(as we have seen in Sect. ||) and M is dominated by the 
evaporation losses. Using Eq. ([!(]) for M. we then find 



A = aQf- 



16V2tt In A ri 



3/2 



Rs 
Ri 



(45) 



where tq is the classical electron radius and 8 — kT c /m c c 2 . 
This result depends on the temperature of the warm disk, 
the viscosity parameter and the distance from the hole. 

We can now make an estimate of the relative extent 
So of the warm disk. Rq is the maximum radius where the 
warm disk can exist (at larger surface density, it devel- 
ops a cool disk component). The optical depth at R is 
thus given by Eq. (^?j) . Computing the optical depth from 
Eq. (O) and equating this to r m we have 



3a c? 



*o 2 



Mi 



1/2 



0.1 ■ M Edd 



(46) 



which shows that Sq is a function of the warm disk tem- 
perature, the viscosity, the relative distance from the hole 
and the accretion rate. Numerically, we have 



So « 0.1 



Rs 



Mi 



M] 



(47) 



Edd 



where we have set InA w 20. With typical values, 8 ~ 
0(1), Ri/Rs ~ 8, Mi/M Edd = 0.1, / = 0.88 and a = 0.3 
this yields So ~ 0.3. We conclude that the extent of the 
warm disk is in an 'interesting' range. It is smaller than the 
distance to the hole, but still large enough that the width 
of the warm disk is large compared to its thickness, so our 
thin-disk treatment of the evaporation region is justified. 

4.4. Conditions for a steady mass flux to exist 

We now estimate the radial distance of the transition ra- 
dius from the hole as a function of the mass accretion rate 
M in the system. At the distance Ri all mass accreted in 
the cool disk is evaporated, M cv = M. The accretion rate 
M of the system is now also equal to the accretion rate 
M in the ISAF and we can write 



M ev w 2 • 2tt RiARi M = 4tt R? S M = M x 



(48) 



We use Eqs. dTo] , ^7| ) and find for Ri (in units of the 
Schwarzschild radius) 



Rj_ 

Rs 



7 x 



,3/2 



_l ( Mi V 

03/2 fl/2 [ MEdd J 



(49) 



Again with the same values as above, 8 ~ O(l), 
M/M Edd = 0.1, / = 0.88 and a = 0.3 this yields 
Ri/Rs ~ 8. We conclude that, using plausible values, the 
estimated transition radius is in the region where we ex- 
pect our model to work. 



5. Discussion and conclusions 

Our model describes a situation which naturally devel- 
ops if a cold standard Shakura-Sunyaev disk is truncated 
within a hot, optically thin flow (ISAF) in the inner re- 
gions of the accreting system. We have proposed a SSD- 
ISAF transition based on few, well-known physical pro- 
cesses: Spitzer's theory of the energy exchange in a fully 
ionized plasma, and standard viscous heating due to fric- 
tion in the accretion disk. We have shown that the tran- 
sition of the cool disk material into an ISAF is the logical 
and inevitable consequence of these basic interactions. The 
process involves two steps: 

(i) At the inner edge of the disk the surface den- 
sity of the cool disk gets low. Virial protons penetrating 
from the ISAF heat the cool disk electrons. The elec- 
tron temperature is limited (T c ~ 80 keV) since they 
can radiate their energy efficiently via bremsstrahlung and 
Comptonization. Once the disk is too thin, proton heat- 
ing overcomes the radiative losses everywhere in the disk. 
The disk heats up, expands and radiative losses become 
even more inefficient. With increasing temperature pair 
production sets in, but also the proton heating gets less 
efficient at higher temperature. Finally a new equilibrium 
for the thin disk is found at several 100 keV. We label this 
the 'warm state', since its temperature is intermediate be- 
tween that of the cool disk and the virial temperature. 
An important aspect in this process is that protons in the 
disk are outside of the main energy channel. The main en- 
ergy is transferred from the external ISAF protons to the 
disk electrons, which loose this energy via radiation. For 
the formation of the warm disk state the internal viscous 
heating of the protons is completely unimportant, but not 
so in the second step of the process: 

(ii) At the temperatures of the warm state the disk 
protons and electrons are not coupled very tightly any 
more, as the timescale for establishing thermal equilib- 
rium is not short compared to the thermal timescale. Now 
the minor energy channel due to viscous heating in the 
warm disk becomes important for the energy budget of 
the protons, because the viscously released energy can not 
be exchanged very efficiently with the ambient electrons. 
The upper part of the warm disk, where the densities are 
lowest and the Coulomb exchange time scales longest, is 
subject to a thermal instability. The size of the unsta- 
ble region depends mainly on the viscosity parameter and 
the temperature of the warm disk. The higher the tem- 
peratures and a, the deeper is the unstable region. The 
protons there are heated to their local virial temperature, 
and become part of the ISAF: the warm state evaporates. 
In the warm state no effective mass condensation takes 
place, since it at the same time becomes optically thin for 
penetrating hot protons. The mass evaporation due to the 
viscous instability of the warm state is therefore the main 
process which determines the mass budget in this region. 

On the basis of this picture, we have also looked at 
the radial structure of an accretion disk in which mass 
exchange with an ISAF takes place. Using the evapora- 
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tion and condensation rates derived, we have investigated 
the conditions under which accretion is possible in such a 
way that the entire accretion flow eventually evaporates. 
In this case a steady state is possible with an inner edge 
to the disk at some finite radius R\ outside the last stable 
orbit. We find that such a steady state is indeed possible 
for plausible values for the accretion rate and viscosity pa- 
rameter. The steady state condition determines a relation 
between the accretion rate and the value of R{. Also, it 
determines the width 5 = AR/ Ri of the warm, evaporat- 
ing disk region; we find J « 0.3 for a of order 0.3 and 
accretion at a tenth of Eddington. 

A potentially important factor which we have not been 
able to include in our picture of the warm disk region is the 
cooling effect of soft photons from the cool disk extending 
just outside the warm disk. If such photons can enter the 
warm region, they will cause a cooling by inverse Compton 
scattering on the warm electrons. Since both the cool and 
warm disk are quite thin (H/R <C 1), the radial optical 
depth of the warm disk is large, and the angle subtended 
by the cool disk as seen from the interior of the warm disk 
is small (to visualize this, see Fig. 5). It thus seems likely 
that the effect of such cooling on the energy balance in 
the warm disk can be small, but this point requires closer 
scrutiny. 

The extent of the region (in distance from the hole) 
where the physical conditions assumed here apply is lim- 
ited. At large radii the proton temperature in the opti- 
cally thin, hot region decreases (oc R^ 1 ). At lower pro- 
ton temperatures the proton penetration depth into the 
cool disk also gets smaller. This limits the region where a 
warm state can be produced. Without the warm state our 
mechanism will probably not work efficiently enough to 
transfer all material from a cool SSD into material from 
an ISAF. Therefore we do not expect our mechanism to 
work at large radii from the black hole. But a combina- 
tion of the coronal evaporation flow model as suggested 
by Meyer et al. (200C| ) and the two stage model proposed 
here could in principle cover a large range in the radial 
direction for the SSD-ISAF transition to occur. 
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